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We propose a scheme for the creation of skyrmions (coreless vortices) in a Bose-Einstein conden- 
sate with hyperfine spin F = 1. In this scheme, four traveling-wave laser beams, with Gaussian or 
Laguerre- Gaussian transverse profiles, induce Raman transitions with an anomalous dependence on 
the laser polarization, thereby generating the optical potential required for producing skyrmions. 
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The recent experimental success in all-optical trap- 
ping of an atomic Bose-Einstein condensate (BEC) |]] 
opens the prospect of studies into the internal struc- 
ture of spinor BECs, including the possibility of creat- 
ing vortex states without core, or skyrmions, in the BEC 
Skyrmions, which do not have an ordinary vortex 
core due to the spin degree of freedom, offer a myriad 
of new physical phenomena beyond those presented by 
other vortex states. One feature is the reduction of ki- 
netic energy associated with the rotation by transferring 
this energy to the spin. Another interesting property of 
skyrmions is that they do not represent a topological ex- 
citation, although their flux is vortex-like far away from 
the line of symmetry j7|. In this paper, we propose an 
optical method to create skyrmions in a condensate of 
alkali atoms. This method, which employs laser beams 
and Raman transition to generate skyrmions in the BEC, 
is related to proposals for the creation of vortices in a 
single-component BEC (|||. 

For a BEC with hyperfine spin F = 1, a skyrmion wave 
function, which is axially symmetric around the z-axis, 
has the form 01 
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with p(z,r±) the total density of atoms, = (x,y) 
the transverse coordinate vector and r± = | r j_ | the ra- 
dial distance from the z-axis. The angle <\> corresponds 
to the orientation of in the o>?/-plane. The func- 
tion (3(r±), which characterizes the spin state of the con- 
densed atoms, is related to the superfluid velocity v s of 
the system by 
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where = (cos (/>, — sin 0, 0) is the unit vector in the (j) 
direction. Because a skyrmion has no vortex core, 



f3(0) = must hold in order to avoid a singularity. For 
/3( r ±) — 7r /2, the superfluid velocity reduces to that of 
an ordinary vortex state. 

To create a skyrmion we consider a BEC which first is 
magnetically trapped in the m = — 1 hyperfine state. We 
assume that the trap is then switched off and the optical 
potential is applied to the BEC. Thus, our initial state is 
given by = yj p(r±, z) and = ^i — 0- O ur objec- 
tive is to design an optical potential which transfers this 
state into the state ([!]). 

Generally a coherent optical potential for atoms is 
created by applying several highly detuned laser beams 
which induce Raman transitions between different an- 
gular momentum states or internal states of the atomic 
BEC (cf Fig. 0). In order to preserve the cylindrical 
symmetry of the skyrmions along the z-axis, the laser 
beams must also propagate along this axis. This restric- 
tion raises the following challenge: the polarization ori- 
entation of laser beams propagating along the z-axis is 
generally presumed to be in a superposition of e x and e y . 
However, as a Raman transition of such beams can only 
effect transfers of the type ip-i to ^i, the possibility of 
creating skyrmions of the form (0) seems to be excluded. 

Fortunately the flexibility in the transverse structure 
of laser beams permits a portion of the laser mode to be 
polarized along the axis of propagation. Mathematically, 
this contribution to longitudinal polarization has its ori- 
gin in the vanishing divergence of the electric field for 
any laser beam; that is, 



divE(r) = . 
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For a laser beam propagating along the z-axis, with 
transverse mode function u(r±) and polarization vector 
e x , the divergence condition (||) makes it necessary that 
the complete positive-frequency part of the electric field 
is given by [[To| 
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where u'(r±) denotes the derivative of u with respect to 
x. The second term on the right-hand side of Eq. (||) is 
essential for inducing transitions to component. 

Our proposal includes four traveling-wave laser beams 
and a homogenous magnetic field B = Be z . The pur- 
pose of the magnetic field is to lift the degeneracy of the 
ground states \m g ) with m g — 0, ±1 being the magnetic 
quantum number of the states. We write the correspond- 
ing Hamiltonian in the form 



H B = E B ^2 m 9\ m g)( m 9 
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The laser beams induce transitions between the F = 1 
manifold of the atomic ground state and an F = 1 man- 
ifold of excited atomic states, such as the D\ line be- 
tween 2 Si /2 (F = 1) and 2 V 1/2 {F = 1) in Rubidium 
( 87 Rb). All four beams are linearly polarized in the re- 
direction and have a large detuning to suppress sponta- 
neous emission. Each beam has a different mode struc- 
ture: we include a Gaussian beam with frequency uoq 
and a Laguerre-Gaussian beam of second order with fre- 
quency U02- Both beams are copropagating along the pos- 
itive z-axis. Our scheme further includes two copropagat- 
ing Laguerre-Gaussian beams (along the negative z-axis) 
of first order with opposite orbital angular momentum 
and frequencies uo\ and respectively. All beams are 
assumed to have the same width w (where w is the width 
parameter for the Gaussian envelope of each beam). The 
electric field of each beam then is of the form ([|), whereby 
the transverse mode functions are approximately given 

by 
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The interaction of the atoms with the laser beams in 
the electric dipole approximation is given by 
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with V is a reduced matrix element. For simplicity we 
assume V to be real. The first contribution describes 
the usual transitions which are induced by x-polarized 
laser beams. The term proportional to u[ has its origin 
from the small z-polarized term which results from the 
vanishing divergence of the electric field. 

The optical potential which can be deduced from such 
a configuration for largely detuned laser beams has the 
general form (see, e.g., Ref. 
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where co res is the resonance frequency for the electronic 
transition. The total Hamiltonian for this system then 
takes the form 

H = Jd 3 x *L 9 <™ 9 |tfi-pK>* m l,+tfni (11) 

m g ,m' g 



with the one-particle contribution 
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where &i is the Rabi-frequency operator for the / th laser 
beam, defined by 

nn ie - iuJit = -d f • E z (+) (8) 
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where m e denotes the magnetic quantum numbers for the 
excited states. Using standard techniques for Clebsch- 
Gordon coefficients (eg, see Ref. ]Tl[), we find for each 
laser beam the expression 
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where = AnTi 2 ai/M is proportional to the correspond- 
ing scattering length a^. Transforming to the interaction 
picture, with respect to Hb of Eq. (||), the optical poten- 
tial assumes the form 
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It is now our task to derive from this general expres- 
sion a special potential that is suitable for the creation 
of skyrmions. As our initial state is the ground- state 
for a BEC trapped in m g = — 1, it becomes obvious 
from Eq. (m) that we require matrix elements of the form 



l|e 2 ^, |0)(— l\e l< P and so on. As the transverse 
mode function u\ is proportional to e*^, its derivative u\ 
is approximately proportional to e 2 ^ -1 ^. This relation 
permits the following construction. 

The detunings = uoi —uj res of the laser beams are the 
most significant frequencies in the problem. To construct 
the desired matrix elements we assume that A2 « Ao and 
A_i w Ai. However, Ao and Ai should be sufficiently 
different so that a Raman transition involving the beams 
and 1 is suppressed. Assuming that the second largest 
frequency is given by the magnetic interaction energy E B 
allows the rotating-wave approximation to be done with 
respect to this frequency in Eq. (14). By assuming that 



5U0Q2 = C^o — ^2 + 2E B <C E B , 00q — Ul>2, (15) 

Soon =uo\— uo-x — E B <C E B , oj\ — lj-i, (16) 
many terms can be neglected, and we arrive at 
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This term can be further simplified because the deriva- 
tive of ui typically scales as 1/w compared to w; itself. 
Including the prefactor 1/k that appears together with 
each derivative, we see that every term with a derivative 
is suppressed by a factor of 
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As the width of a laser beam always exceeds its wave- 
length A, the factor ( |l8| ) is always much smaller than 
one. For instance, for an optical transition (A = 795 nm) 
a quite strongly focused beam with a width of 5 jam has 
e « 0.025. This allows us to develop V op t in a power 
series in e. Before doing so we observe that in Eq. ( |l7| ) 
the transitions between |0) and | =L 1) are suppressed by 
a factor of e as compared to the transition between |1) 
and I — 1). Since this is undesirable for the production of 
skyrmions we assume that the intensity Aq (and hence 
uo) of the Gaussian laser beam defined by Eq. ^ will 
be of the order of e as compared to the other prefactors 



A\. Inserting the mode functions (^|) into Eq. ( |l7| ) and 
developing Vo P t to first order in e, we find 
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Eq. ( p^9|) is the main result of this paper. It consists of 
several parts which will now be analyzed. 

The first row in Eq. ( |l9l) represents an optical trap- 
ping potential if the laser detuning Aj is positive. For 
f = r±/w <C 1 it is harmonic. Its physical origin is 
a Raman transition which returns to the same internal 
hyperfine level. As the state \m g = 0) has two possibili- 
ties (\m e = ±1)) to make such a transition the trapping 
potential is twice as strong as compared two the states 
\m g = ±1). Note also that this term is of order e° and 
therefore represents the dominant contribution. 

The second row describes ordinary Raman transitions 
between the hyperfine states \m g = — 1) and \m g = 1). 
As laser 2 corresponds to a Laguerre-Gaussian beam of 
second order, an orbital angular momentum of 2% will 
be transferred to the atoms, as it is required for the for- 
mation of the skyrmion (|l|). We remark that the same 
result could have been achieved by replacing the Gaus- 
sian beam by a weak Laguerre-Gaussian beam of order 
1, and beam 2 by a Laguerre-Gaussian beam of order 
— 1. The optical potential would only be changed in that 
the trapping potential in the first row of Eq. (|l|) would 
become harmonic (apart from the overall exponential fac- 
tor). From an experimental point of view this alternative 
may be advantageous because it only requires two differ- 
ent types of laser beams (Laguerre-Gaussian of order 1 
and -1) instead of four types. 

The third and the fourth row in Eq. ( fL9l) describe 
anomalous Raman transitions in which the atoms ab- 
sorb or emit a photon with polarization e z in one of 
the two processes involved. In Fig. |l] this corresponds 
to a vertical transition from \m g = 1) (or \m g = —1)) 
to \m e = 1) (or \m e = —1)). The possibility of these 
transitions is a direct consequence of the full mode struc- 
ture (||). Since these vertical transitions are proportional 
to the derivative u[ of the transverse mode functions, 
no orbital angular momentum is transferred even though 
the corresponding beams are of Laguerre-Gaussian type 
±1 (more precisely, for r <C 1 the terms proportional 
to ex.p[±2i(j)] in Eq. ( [l9| ) are negligible). The second pro- 
cess which completes the Raman transition is an ordinary 
transition from \m g ) to \m e = m g ±1). Since the cor- 
responding laser beam is of Laguerre-Gaussian type ±1 
the total orbital angular momentum transferred to the 
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atoms is ±7i. As a result, the anomalous Raman transi- 
tions create simple vortex states in \m g = 0) and doubly 
excited vortices in \m g = 1). 

It might be counterintuitive that the z-polarized con- 
tribution to the electric field (||), which usually can safely 
be neglected, does now play an important role for the cre- 
ation of skyrmions. In fact, its contribution remains as 
small as usual compared to the dominant parts (diagonal 
elements of V opt ) . The important point in our proposal is 
that we can make all matrix elements inducing internal 
transitions equally large by using a weak laser beam so 
that the usual Raman transition rn g = — 1 <-> rn g = 1 is 
not stronger than the anomalous Raman transitions. 

As the diagonal elements of V opt are much larger than 
the off-diagonal elements they provide a strong trapping 
potential. According to Eq. (|l|) a skyrmion is a super- 
position of the ground-state in m g = — 1 and a singly 
(doubly) excited vortex state in m g = (m g = 1), re- 
spectively. Because of the strong trapping potential these 
states all have different energies. Thus, to induce transi- 
tions using the off-diagonal elements of V opt these matrix 
elements need to rotate at some frequency to match a 
resonance condition. This can be done by adjusting the 
frequencies Slj 02 and Suju appropriately. The efficiency 
of the skyrmion production strongly depends on 5wo2 and 
Suju. 

We have numerically examined the time evolution of 
a spin-1 BEC in the potential ( |19|) . Our simulations are 
based on the split-step method Jl4| in two spatial dimen- 
sions for a three-component nonlinear Schrodinger equa- 
tion (It is assumed that the BEC is homogeneously ex- 
tended along the z-axis with a length of 50 /im). We con- 
sider a BEC of 15000 Rubidium atoms (M = 1.45 x 1(T 25 
kg, a s = 5.4 nm, a a = —0.05 nm Q). For a much larger 
number of atoms our method seems to be unsuitable be- 
cause the potential (|l9| ) has a maximum height due to 
the exponential envelope. A strong repulsive interaction 
between the atoms then pushes the atoms out of the trap- 
ping area. 

For the optical potential (|l9| ) we consider the situation 
that Al/Ao = 2Al/Ai holds and the amplitude of the 
Gaussian beam is given by Aq — 3eA2 / y/2. The effective 
Rabi frequency of the trapping potential is then given by 
V 2 AI/(ATi 2 Aq) which we take to be 66 kHz (For a de- 
tuning Ao of 1 GHz this corresponds to a laser intensity 
of about 25 mW). The width of the laser beams is 5 /im. 
For the results shown in Fig. [|we have used 5ujq2 = — 2.4 
kHz and 5uo\\ = 2.95 kHz. At time t = 6.6 ms this re- 
sults in a superposition of internal states with about 27% 



of the atoms in state m a 



-1, about 34% in the state 



m g = 0, and 39% in state m g = 1. For various other 
choices of Suoo2 and Suju we have found as little as 4% 
population in the initial state m g = — 1 and as much as 
70% in one of the other two states. The phase of the three 
components agrees very well with the situation described 
by the state (Q). Animations of the time evolution are 
available on the internet |H| . 
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m g= -l 



FIG. 1. A scheme of Raman transitions needed for the creation of skyrmions. The solid line corresponds to laser beam 0, the 
dashed line to beam 1, the dot-dashed line to beam -1, and the dotted line to beam 2. Diagonal transitions describe ordinary 
transitions with photons polarized in x-direction. Vertical transitions are anomalous transitions with photons polarized in beam 
direction (2- axis). 
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FIG. 2. Time evolution of a three-component BEC under the influence 
the three hyperfine states m g — 0, ±1 at four instants of time (given in ms) 



The initial state is the ground state (of a harmonic 



trap) in the state m g 



-1. 
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